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A theory of flame propagation in curved channels is developed within the framework of the on- 
shell description of premixed flames. Employing the Green function appropriate to the given channel 
geometry, an implicit integral representation for the burnt gas velocity is constructed. It is then used 
to derive an explicit expression for rotational component of the gas velocity near the flame front by 
successive separation of irrotational contributions. We prove that this separation can be performed 
in a way consistent with boundary conditions at the channel walls. As a result, the unknown 
irrotational component can be projected out by applying a dispersion relation, thus leading to a 
closed system of equations for the on-shell fresh gas velocity and the flame front position. These 
equations show that in addition to the usual nonlocality associated with potential flows, vorticity 
produced by a curved flame leads to specific nonlocal spatial and temporal influence of the channel 
geometry on the flame evolution. To elucidate this influence, three special cases are considered 
in more detail - steady flame stabilized by incoming flow in a bottle-shaped channel, quasi-steady 
flame, and unsteady flame with small gas expansion propagating in a channel with slowly varying 
width. In the latter case, analytical solutions of the derived equations are obtained in the first 
post-Sivashinsky approximation using the method of pole decomposition. 

i 

I. INTRODUCTION 

CO 

o _ 

Analytical description of premixed flame propagation is one of the most difficult and exiting tasks in combustion 
^•-y science. In its full generality, it represents a tremendously complicated problem involving several length scales, nonlocal 
[ interactions, and various instabilities inherent in the governing mechanisms, which in turn necessitate consideration 
Oh. of flame dynamics in the regime of saturated nonlinearity. On the other hand, this problem naturally admits several 
important simplifications which make it amenable to theoretical treatment. First of all, the fact that the deflagration 
— , , is an essentially subsonic process allows the gas flows be considered incompressible, both in the fresh and burnt 
' gas regions. Also, neglecting viscous effects in these regions, gases can be treated as ideal outside the flame front. 
These assumptions turn out to be very good approximations to real flames. Furthermore, large separation of scales 
determining inner dynamics of the flame front and its global behavior allows the use of asymptotic methods to derive 
local relations between the fresh and burnt gas flows at the flame front. These relations contain information about 
transport processes inside the front, and are expressed in the form of jump conditions for the gas velocity and its 
pressure across the front, and an equation determining local consumption rate (the evolution equation). They play 
the role of boundary conditions for the large-scale flow at the front considered as a surface of discontinuity. 

Yet, even simplified this way the problem is still exceedingly complicated. It requires solving a system of nonlinear 
partial differential equations on both sides of the flame front, to be chosen so as to satisfy the jump conditions across 
the moving front. It is clearly impossible to resolve this system explicitly, because even if the flame propagates in an 
k>( \ initially quiescent fluid, so that the gas flow is potential upstream, this is no longer true downstream the flame because 
of vorticity generated by the curved flame front. In this connection, it is important to note that flames used in practice 
are characterized by large thermal expansion of gases: the fresh to burnt gas density ratio, 9, is normally 5 10. It is 
not difficult to show [l| that under such circumstances, the weak nonlinearity approximation is not applicable to the 
description of flame evolution (except the earlier stage of development of the Darrieus-Landau instability @, Q), so 
that the problem of solving the flow equations is faced in its full generality. 

This difficulty precludes theoretical analysis if one is concerned with the explicit structure of the burnt gas flow. 
However, it is evolution of the flame front, its position and shape, which is of greatest interest in practice. As was 
shown in 0, @, 0] > this limitation of the problem opens a way round the above-mentioned difficulty. Namely, it 
turns out to be possible to derive an exact system of equations describing the flame-front evolution, which is closed 
in the sense that it involves only restrictions of the flow variables to the front (their on-shell values). Not only this 
approach permits the well-known results on the flame-front dynamics be derived anew in a simple and elegant way, 
it also brings to the scope of theoretical treatment classical problems which have not been accessible by means of the 
conventional analysis [a, H- 

Considerations of [IlUJa S B El have been restricted to the simplest case of flames propagating in straight channels. 
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A question of principle which is very important from both experimental and theoretical points of view is whether the 
results obtained can be extended to channels with curved walls. Importance of this question lies in the fact already 
mentioned that the flame evolution is an essentially nonlocal process. In view of this nonlocality, flame propagation 
is naturally expected to be affected by the wall curvature. This effect is anticipated to be especially pronounced in 
the two-dimensional case due to peculiar long-range behavior of the Green functions. 

The purpose of this paper is to show that in the two-dimensional case, the above question resolves in the positive. 
Namely, the on-shell flame description will be extended to symmetric channels with arbitrarily curved walls. The 
influence of channel geometry on flame dynamics turns out to be much more complicated than predicted by the 
potential flow models [l(| El]. Because of vorticity produced by the curved flame, this influence does not reduce 
to a mere conformal transformation of the flow variables. Apart from this fact which is clear from the outset, our 
investigation reveals specific spatial and temporal nonlocalities associated with the variable channel shape. 

The paper is organized as follows. Derivation of the main integro-differential equation for the on-shell gas velocity, 
called for brevity the master equation, proceeds along the lines of Refs. 0, Q. We first derive in Sec. HJan implicit 
integral representation for rotational component of the burnt gas velocity using the Green function approach, and 
discuss its validity for various asymptotic channel geometries. In Sec. 1111 A| we generalize the dispersion relation for 
irrotational velocity component to arbitrarily curved channels. This consideration identifies conditions under which 
irrotational velocity can be safely neglected, which is then used in Sec. IIII Bl to gradually simplify the integral formula 
for the rotational velocity component. An explicit expression for this component at the flame front is obtained in 
Sec. IIII C[ which is the last step of the derivation of the master equation written down in Sec. IIII Dl To elucidate the 
structure of this equation, we consider three special cases - steady flame stabilized by incoming flow in a bottle-shaped 
channel (Sec. lIVp . a quasi-steady flame, and a flame with (6 — 1) <C 1 propagating in a channel with slowly varying 
width (Sec. EJ, which admits complete analytical investigation. The results of the work are discussed in Sec. I VII 

II. INTEGRAL REPRESENTATION OF THE FLOW VELOCITY 

Consider a 2D flame propagating in a symmetric channel filled with an initially quiescent uniform ideal reactive 
fluid. The channel width is assumed to vary smoothly along the channel, but is otherwise arbitrary. Let the Cartesian 
coordinates (x,y) be chosen so that the y-axis is along the symmetry axis, y = — oo being in the fresh gas. Then the 
flame-front position at time t can be described by an equation y — f(x, t). The fluid velocity (w, u) will be measured 
in units of the planar flame speed relative to the fuel, and the fluid density in units on the fuel density, 9 > 1 denoting 
its ratio to that of the burnt gas. In the subsequent analysis, we shall widely use complex notations combining the 
coordinates and velocity components into the complex variables z = x + iy and u> = u + iw. In addition to the physical 
z-plane, it will be convenient to introduce an auxiliary £-plane, C = ? 7 + with 77, £ real, such that the function 
£ = g(z) maps the curved channel in the physical plane onto a straight channel of constant width, 25, in the auxiliary 
plane. In coordinate form, this mapping will be written simply as rj — rj{x, y), £ = £(x, y), and similarly for its inverse. 
We assume that g{z) is analytic on some open subset in z-plane, containing the physical channel, and that g(iy) = i£, 
£ = — 00 corresponding to the fresh gas (see Fig.[T|). The symmetry of the physical channel with respect to the y-axis 
implies that the function g(z) possesses the following reflection property 



g(-z) = -g(z), (1) 

where the bar denotes complex conjugation. Accordingly, we will be concerned with symmetrical flame configurations 
satisfying 

f(x,t) = f(-x,t), w(x,y,t) = -w(-x,y,t), u(x,y,t) = u(-x,y,t). (2) 
Under the assumptions that the gas flow is ideal incompressible, its velocity obeys the following equations in the bulk 



0, (3) 



dw du 
dx dy 

da da da ftS 
!H +W d-x +U d-y= > (4) 



where 

rr = 

dx dy 



du dw 

a = ( 5 ) 
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is the vorticity. Let us introduce the stream function ip according to 

dtp dip ,^ 

Then the continuity equation §S§ becomes an identity, while the definition ([5]) takes the form 

d 2 d 2 

Aip = a , A.- + -. (7) 

In order to construct an implicit integral representation for the burnt gas velocity, let us treat the latter relation as 
the Poisson equation for tp, considering the vorticity distribution downstream of the flame front as a given function. 
Since the stream function must be constant along the channel walls, we have to find the Green function, G(x, y; x, y), 
of the Laplace operator, satisfying 

AG(x, y; x, y) = S(x - x)S(y - y) , G(x(±b, £), y(±b, £);£, y) = const , 

where S(-) is the Dirac function. It is not difficult to verify that one can choose 

G(x,y;x,y) = —hi tfjL- _J + c .c., (8) 

'{ S [ff(*)+*(*)J} 



where "c.c." stands for the complex conjugate of the preceding expression. Indeed, the right hand side of Eq. JHJ) is 
analytic everywhere in the channel except the point z = z, therefore, it satisfies the Laplace equation. On the other 
hand, in a vicinity of z = z 

G(x, y; x, y) = — ln(z — z) + R(z, z) + c.c, 

47T 



where R(z, z') is analytic at z = z' . Hence, denoting r = y/Jjx — x) 2 + (y — y) 2 = \l [z — z)(z — z), one has 

G(x,y;x,y) = ^rlnr plus terms satisfying the Laplace equation, so that G correctly reproduces the ^-singularity 
at (x,y) — (x,y). Finally, it is not difficult to check that G(x(±b, £), y(±6, £); x, y) = 0. Thus, the general solution of 
Eq. |(7|) can be written as 

i 1 (x,y)= dsG(x,y;x,y)a(x,y)+ipo(x,y), (9) 
is 

where S denotes the (open) region downstream the flame front, ds = dxdy, and %po(x,y) the general solution of 
Aip = 0, satisfying the boundary conditions 

tpo(x(-b,^),y(-b,^)) = const_, ip (x(+b 7 £), y(+b, ^)) = const+ . (10) 

The difference of the constants const± gives the overall mass flux of the burnt gas. Following the usual procedure 
013 1 we now rewrite the first term in Eq. ([9]) as an integral over fluid particle trajectories. Namely, at each given 
time instant t, we consider E as spanned by fluid elements that crossed the flame front during its evolution up to the 
given t. Symbolizing differentiation with respect to x by a prime, let us introduce the quantities N = (1 + J' 2 ) 1 / 2 and 
= (v + ,n), where 

v + = (w + ,u + ), u+(x,t) =u+(x,t) — — , 

n is the unit vector normal to the flame front (pointing towards the burnt gas). The subscript "+" ("— ") is used 
throughout to indicate restriction to the front of the corresponding downstream (upstream) flow variable; thus, w" 
is the normal burnt gas velocity relative to the flame front. Then, taking into account that the "volume" ds of an 
element is conserved because of flow incompressibility, as is the vorticity [Cf. Eq. Q], we can write 

dsG(x,y;x,y)a(x,y) =tp v (x,y) (11) 

is 

r* 1 mn{£r\g(z)-g(Z{x,t,T))]} 
dx / drN(x, t)v$ (5, T)cr+(i,r)— In rV + c.c. 



B J-oo 47T cog , 



{5 [9(z)+9{Z(x,t,T))}} 
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where B — B(t) is the abscissa of the right end-point of the flame front (it is the solution of r](B + if(B,t)) = b), 
and Z(x,t,T) = X(x,t,r) + iY(x,t,r) the current position (in z-plane) of a fluid clement that crossed the point 
(x, f(x, t)) at the flame front at time r <t. Using Eq. fl}, and taking into account that <r+, X are odd functions of x, 
while N,i>™, and Y are even, it is not difficult to see that ip v (—x,y) — —^) v {x, y), so that the corresponding velocity 
field satisfies Eq. ([2]), and therefore, so does the irrotational contribution generated by tpo. 

To ensure convergence of the integral over r, one has to restrict the class of allowed functions g(z), and/or introduce 
some assumption about the flow at infinity. The simplest choice is to assume that the channel is asymptotically 
straight, i.e., g{z) — > const ■ z for Imz — > +00, with an initially planar flame propagating from infinity. Then ct + (x,t) 
is exponentially small for r — > —00. However, this assumption is very strong and can be relaxed. In fact, it may be 
sufficient to have an asymptotically V-shaped channel, with the opening angle 2a < n. Indeed, in this case 

1 \ ib . z 
g{z) = — in — , 1m z — > +oo, 
a i 

so that for large |r|, the logarithm in the integrand of Eq. (TTTj) is proportional to 

-- MtiooM) - 2ig(z)] } ~ Ko|t|)-^ 2q , 
a 



ex ^Tb 



where u^, — dY(x,t,r) jdt. Thus, if a < ir/2, the r-integral converges (since the other factors in the integrand are 
bounded functions of t), provided that Uoc does not vanish for the gas elements carrying nonzero vorticity. But the 
latter condition is just what is normally observed in experiments with Bunsen flames, which indicate that vorticity 
produced by the flame is carried away in thin jets spanned by the stream lines with ^ 0. Restrictions as to the 
asymptotic channel geometry can be weaken even further taking into account parity properties of the flow variables. 

The last comment concerns analytical properties of the flow variables. As always, we assume that the flame is stable 
with respect to short wavelength disturbances i.e., that there is a short wavelength cutoff, A c , whose precise value 
is determined by specific transport processes inside the flame front. This cutoff ensures smoothness of the functions 
involved, which are considered throughout as infinitely differentiable. 



III. DERIVATION OF THE MASTER EQUATION 

A. Dispersion relation for irrotational velocity component 

It is impossible to compute the right hand side of Eq. (jTT|) in general, because for that one would had to solve the 
system of hydrodynamic equations downstream explicitly. However, in order to describe dynamics of the flame it is 
sufficient to determine the flow structure in a vicinity of the front only. This can be done by decomposing the burnt 
gas velocity as 



where v p is a potential component incorporating the unknown information about the bulk flow, which is eventually 
"projected out" using a dispersion relation. To derive this relation, we require v p to satisfy the following conditions: 

a) divw p = 0, iotv p = 0, implying that uj v = u p + iw p is analytic function of z € S, 

b) the function fi(£) = u) p {g- 1 {Q) d9 ^ is bounded in #(£), 

c) fi(— b + i£) = fl(+b + i£) for all £ downstream. 

Consider the region downstream of the flame front in the auxiliary plane. Since uj p {z) is assumed analytic downstream, 
0(C) is analytic for £ G so that one has, by virtue of the Cauchy theorem, 

i/ as(s) ^ COt {| (< ~- C) } 0(C " ) = fi(C) ' (12) 

where <9g(£) is the boundary of g(£), ran counterclockwise by £. If 0(C) satisfies the condition c), the wall contributions 
to the contour integral cancel each other because of the 2&-periodicity of the cotangent. If, in addition to that, b) is 
also met, then taking the limit ( — > C+ = 9( z +)> z + = x + if(x, t), in Eq. (jT2J) gives 



If „ r 7T „ -1 1 f+b+ioo _ _ 

X dC cot {-c-c+l o(c) -— / dcfi(c=n + , 

2bi J g(F) 1 2b J 2bi J_ b+ioc 
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where the integral over the image g(F) of the flame front is understood in the principal value sense. We can get rid 
off the constant contribution of the infinitely remote interval by differentiating this equation with respect to x. After 
that, integration of the first term by parts taking into account condition c) yields 

or equivalently, 

fetl j + * dxco t[^[g{2) - g( Z+ )]) (fi + )' = (fi+)' , (14) 
where z — x + if{x, t) , fi + = £l(g(z)). 



B. General structure of rotational velocity component 



We now proceed to the derivation of an explicit expression for the rotational velocity component, v v . Let the 
equality of two functions tpi(x,y) : tp2(x,y) up to a field satisfying a) - c) be symbolized as <pi = <p2- Then Eqs. ([5]), 
(fTTj) give for the complex burnt gas velocity 



4tt 



tlx 



f 7r 

+oo sm\ — {g(z) - g(Z(x,t,t 

d T M(x,t-T)\n- L4 ° 



"))]} 



; {5 [9i z )+9(Z(x,t,t-T))]} 



+ c.c. 



(15) 



where we denoted $ = d/dx — id/dy, introduced the memory kernel M(x,t) = N(x,t)v™(x,t)a+(x,t), and changed 
the integration variable r — » t — r. In fact, condition a) is obviously fulfilled by the velocity vq corresponding to the 
term ipo in Eq. @. Furthermore, under the assumption about convergence of the r-integral, discussed at the end 
of Sec. m the right hand side of Eq. (|15p is bounded in E, and therefore, so is the velocity ujq = uq + iwo, because 



the total burnt gas velocity is bounded. Hence, the function O (C) = u (g 1 (Q) 
condition b) is met. Finally, to show that c) is also fulfilled, we write 

d 5 _1 (C) _ -dx dy 



dg-HO ia 



is bounded in g(E), so that 



and then 



Q = u 



dy 



dx 



w - 



dy_ 



dx 



u - 



Differentiating the boundary condition (| 1 0[) with respect to £, one finds 

w o|r ~ "0^7 J (±6, = , 



'0£ 

which shows that the imaginary part of Slo vanishes at the walls. On the other hand, differentiation of the reflection 
property (fT]) written in the form x(— rj, £) = — &/(— f?, £) — y(r],^), gives 

Together with ito(— x,j/) = uo(x,y), Wo(—x,y) — —Wo(x,y), these relations prove that Qo(— & + = flo(+b + 
Next, consider the kernel of the i-integral in Eq. (TTf 



K(x,t,Q = ± 



+00 sin 
drM(x,t- t) In 



{^{t-g(Z(x,t,t-T))}} 



cos • 



{lb [< + 9(Z(x,t,t-T))]} 



c.c. 



(16) 



Premixed Flame Propagation in Curved Channels, v 1 - October 13, 2009 



6 



Our task is to evaluate this function in the case when the observation point z is close to the flame front. To that end, 
we have to perform first some general transformations on K, which is the subject of the present section. Before we 

proceed to the computation, let us prove the following simple but important 

Lemma: if the complex- valued function p = p(x, /i) (/i is a parameter) satisfying p(—x, (m) — —p(x, fj,), is such that the 
point g(p) £ <?(£), and if m(x) is odd, then the field 



B sin\ — [g{z) - g(p)]\ 

dxm(x) In ■ 1 46 J 



{J[fl(*)+5(P)]} 



+ c.c. 



(17) 



satisfies conditions a)-c). Indeed, noting that for ip(z) analytic, 

d(p(z) 



dz 



, v(«) = o > 



(18) 



we see that uj p (z) is analytic in S. Furthermore, the corresponding function fl p — tu p (g 1 (£))dg 1 /dC, takes the form 

d 



n p (C,«i,r) 



d( 



- B sin l(-g(p)}\ 

dxm{x) In ■ 1 4 " ' 



C0S {Tb[ c + g{p) ]} 



It is bounded for £ S <?(£), so that b) is also met. Finally, c) is easily verified changing the integration variable 
x —x, and taking into account oddness of the function m(x). As a consequence of this result, differentiation of 
Eq. (TrT|) with respect to /i shows that the function 

' /7 «m(«) { l K - , Wl} ^ - tan {J [< + Wj] } ^) 

also satisfies a)-c). 

We now rewrite the kernel integrating it by parts. For this purpose, we introduce the quantity 



/ 



dnM(x, t-ri)= M(x, t, t) 



(19) 



and abbreviate Z(x, t,t — r) = Z T , g(Z T ) = g T . This yields 



K(x,t,0 = M(x,t,0)— In 



sin{J[ C - 9 (f)]} 



C0S {^ [ C + ff(5) ]} 
drM ( i ' t ' r ) ( cot {^ [C-.9r]}.9r - tan [C + g T }}t T ) + c - c - 



(20) 



where the dot denotes differentiation with respect to r. By virtue of the Lemma, the first term here can be omitted, 1 
because E is open, while z £ F<f2,. 

Let us allow r to take complex values. Under the assumption of existence of a short wavelength cutoff, all flow 
variables are smooth functions of time. Hence, these functions (in particular, the memory kernel M) are analytic in a 
vicinity of the real axis in the complex time-plane. Therefore, the only singularity of the integrand in Eq. (|20j) . which 
can cross the real axis, is at r such that 



g{z T ) = C- 

Specifically, the imaginary part of r vanishes at the point {r, x} satisfying 

X[x, t, t — t] = x , Y[x,t, t — t] = y . 



(21) 
(22) 



1 This will be symbolized by the same sign = used previously to relate velocity fields. 
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The pole corresponding to r satisfying this equation determines the structure of the rotational component of the 
burnt gas velocity. We extract the pole contribution by deforming the contour of integration over r. Specifically, we 
write J dr = 1/2 J dr + 1/2 J dr, and deform the contour of integration in the two halves symmetrically with respect 
to the real axis, moving it pass the solution of Eq. (|2"2"|) , for all x (see Fig. [5J where C+ and C_ denote the upper 
and lower branches of the contour, respectively). Since the integrand in Eq. (|20[) involves g T , in order to promote r 
to a complex variable the operation of complex conjugation is to be further specified. It is convenient to define g T 
as an analytic continuation of (g T ) with respect to r (in other words, the argument r in ~g T is not subjected to the 
conjugation). The pole contribution to the integral kernel then reads, by virtue of the Cauchy theorem, 

i 

K (x,t,C) = -x(Imr )M(i,i,T ), 

where x( x ) is the sign function, 

ami to = To(x,t, f) denotes the solution of Eq. (|2"Tj) . To keep the boundary conditions satisfied (see condition c)), we 
must perform the same contour deformation for both terms in the integrand in Eq. (|20|) . But when the observation 
point is at the channel wall (say, the right), one has g T0 = +b + if, hence £ + g^- — 2b, which is a pole of the tangent. 
In order to guarantee continuity of the extracted contribution, we must proceed with contour deformation so as to 
embrace the pole of the tangent, again for all x and (eE. Let us denote this pole ctq = ao(x,t, f). Thus, 




a (x,t,b + i£)=To(x,t,b + iQ. (23) 

Furthermore, the change x — > —x transforms Z T — > —Z T (r is not conjugated), hence, using the property |T]), 
g T — > — g T , which together with Eq. (|21|) imply 



r (-£,*,-& + if) = To(x,t,b + i£). 

Combining this with Eq. (f2"3"]l gives 

t (-x, t, -b + if) = a (x, t,b+ if), (24) 

In other words, inversion x — > —x maps the ro-poles at the left wall to the <7o-poles at the right wall. Let us denote 
<7_i = <t_i(x, t, f) the pole of the tangent which is in the same correspondence with the right-wall value of To, viz., 

T (-x,t,b + iQ = a-i(x,t,-b + iO- (25) 

Again, continuity requires the contour to encompass the a_i-pole. In turn, this pole merges at the right wall a pole 
of the cotangent, to be denoted r_i(i, t, f), with the correspondence given by the inversion 

r_i -6 + if) =o-_! (£,*,& + if). (26) 

Proceeding in this way, we find that all poles of the integrand in Eq. PD)) must be embraced by the contour, the pole 
enumeration being defined by the following infinite chain of relations for their boundary values 

T n (-x,t,-b + i£) = <r n (x,t,b + i£), 
T n (-x,t,b + i£) = a n -i(x,t,-b + i£), (27) 

for all integer n. Note that for all x, each a- or r-pole, except To, belongs either to the upper or lower half of the r-plane. 
In fact, vanishing of the imaginary part of any of these poles for some x would mean that there is a real trajectory 
(X,Y) connecting the point (x,f(x,r)) with the observation point (x,y), which is impossible since (x,y) S S, while 
all the image sources corresponding to the r- or a poles are outside of the channel. It is not difficult to see that Im tq 
is negative (positive) when the observation point is at the right (left) wall. It follows then from Eqs. (f2"4"]) - (|2"T|) that 

Imr ra = Im<7„ = x( n )> n 0; Imao > 0. (28) 

Thus, the integral kernel takes the form 

+oo 

K(x,t,C) = |x(ImTo)M(f,t,ro) - ^M(x,t,a ) + % - ^ x(n) {M(x,t,r n ) -M(x,t,a n )} 

n— — OG 

-TQbJ cuc drM(x,t,T) (cot{ J[(- 5T ]}5 T -tan{^[C + g r ]}5 r ) + c.c. (29) 
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Integration over r <E (— oo,0) in the last term of this expression gives rise to a contribution which is = 0. Indeed, 
these t's correspond to the integration over trajectories before they crossed the flame front (it is meant that these are 
trajectories of the burnt gas elements, continued to r < 0, not fresh gas elements). In other words, Z T ^ E for r < 0, 
and the Lemma allows us to omit this contribution. Furthermore, for all r's corresponding to the remote semicircles, 
excluding a finite domain near some of the remote r- or er-poles, the arguments of the tangent and cotangent have large 
imaginary parts of the same sign, hence, cot(-) — > ±i, tan(-) —> ^i, and the integrand becomes ±iM(5, t, r) (g T + g T ), 
plus an exponentially small remainder. The leading term is canceled by its counterpart from "c.c." Indeed, using the 
definition of the function g T , we find 

iM{x,t,T)(g T +t T )+c.c. = iM{x,t,r) ( ^ + ^] - iM{x,t,r) 1 dgT 



dr dr J \ dr dr 

Changing the integration variable r — > r in the second term, and taking into account that C + , C_ are complex 
conjugates of each other shows that this expression is zero. Thus, the integral over the expanding semicircles vanishes 
exponentially. On the other hand, if the sum in the first line of Eq. (|29[) converges, the contribution of the excluded 
domain also vanishes. In this case, therefore, the integral in Eq. (f2"9"| can be completely omitted. 

In the course of the contour deformation, one may also encounter singularities of the functions M(x,t, r), g T 
themselves. Under quite general assumptions about their analytic properties, it can be proved that these singularities 
contribute terms = 0. Below we consider only the case when M(x,t — r) is a meromorphic function of t, i.e., it is 
allowed to have any number of poles of arbitrary order. Any /cth-order pole of the function M(x,t — r) becomes 
a (k — l)th-order pole of M(x, r, i). Consider the case when M(x,t — t) has a second-order pole at some complex 
t = t* (t* is generally a function of x and t, but for brevity, we omit these arguments; also for definitcness, we assume 
that ImT, > 0). Then M(J,t, t) contains a term m{x,t) / '(r — t„), with some m(x,t) satisfying m(— x,t) = —m(x,t). 
Reality of the function M(i, r, t) implies that it has also a pole at r = r,, and that the two poles have the same real 
residue. Crossing the pole r* by the contour C+ adds to the right hand side of Eq. (|29|) a term 

2m 
^ 165 



m{x,t) (cot [C-StJ}st, -tanj^j [C + g T J } 5 r .) + c.c. (30) 



We see that this brings in an infinite sequence of singularities to the kernel K(x,t,Q considered as a function of £. 
Namely, when £ is close to g Tt + 4bl, I € Z, the new term reduces to a simple pole 

im(x,t) g Tt 



2 C - 9t, - 46Z ' 

On the other hand, since the poles of the cotangent are in one-to-one correspondence with integer numbers, there is n 
such that g Tn = £ — Abl. It follows from Q — > g Tr +4bl that r n — > r*. In particular, Imr„ > 0, therefore, the nth term of 
the series on the right of Eq. ([29]) contributes +(i/2)m(x, t)/(r n — t*). Furthermore, g r?l = 5r» +ffr» (^n — t*)+o(t„ — r*), 
or, g rai = C — 46Z — g Ti (r„ — t») + o(r„ — r*), so the sum of the two contributions is 

im(x, t) g Tt i m(x,t) = 

2 C~9r,~ 4bl 2t„-t, U ' 

A similar consideration shows that the same is true of the pole r* , and of the a-contributions (the latter cancel the 
poles of tangent in the expression (f30|) ). When I runs over all integers, so does n, and we see that all pole contributions 
to the right hand side of Eq. (|2"9"|) cancel. Differentiating with respect to r*, it is not difficult to see that poles with 
k 3 do not contribute either. 

Next, let M(x,t — t) have simple poles at r*, r*. These poles contribute to the function M(x, t, t) a term 
m(x,t)ln{(r — t*)(t — r*)} , with some m(x,t) having the same properties as before. Crossing the logarithmic sin- 
gularities leads to the 2ir jump in arg(ln(-)) for all points of the contour C_ U C + , located at one side of the point t* 
or r* , so the pole contribution to the integral kernel reads 

^m(x,t) [ dr (cot { ^ [C ~ 9r\ } 9t - tan j ^ [( + g T ] | 7^ r ) + c.c. 

= --m(i,i) [lnsinj— [C-.g r ]| -lncosj— [C + ff r ]|J_ + c.c. (31) 

(we may assume that the contour avoids singularities of the integrand, because they are isolated. Then all singularities 
of the integral are associated with the contour endpoints). For £ close to g Tt + 46/, I G Z, the singularity is 

im(x, t) , „ 
V^ ln (C-.9r, -460- 
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It is compensated by +i/21n(r„ — r*) coming from the n-th term of the series in Eq. (|29)l . where n is such that 
g Tri = ( — 46L Again, cancelation of singularities takes place for all n. 

Thus, we arrive at the conclusion that the contour deformation "wipes out" all poles of M, so that no pole or 
logarithmic singularity is left in the kernel K. 

Now, this result can be used to prove that any pole or branch singularity in g T gives rise only to terms = 0. Namely, 
if g T has a pole at r* , then its contribution to the right hand side of Eq. (|29[) reads 



M.(x, t, t* 
166 



j dr (cot j-j^ [C-5t]}st -tanj^ [C + 9 T ]jg T ) +c.c, 



(32) 



where C* is a contour surrounding t», which can be taken as small as desired. Furthermore, if there is a finite cut in 
the r-plane, connecting branch points of g T , it brings in a term 

~16& J c drM ^' t ' T ^ ( cot (i^ iC-9r]^9r -tanj^ [C+g T ]jtr) + C -C-, (33) 

where Co is a contour embracing the cut. Both expressions (|32p and (|33p contain logarithmic singularities with respect 
to £. But the right hand side of Eq. (|2"9"|) has just been proved to be free of such contributions, so the new singularities 
have no counterterms. At the same time, the burnt gas velocity is bounded, as are all the irrotational contributions 
omitted in the course of its simplification [Cf. condition b)] . Hence, the singular point and the cut must be such that 
g T ^ g(E) for t = t* or t £ Co. Noting also that since M is odd in x, so is [M], and applying the Lemma, we see that 
the above expressions = indeed. 

To summarize, under the assumption of convergence of the series over r- and cr-poles, the function (|16p can be 
written 



K(x,t,0 =K c (x,t,0 + K c (x,t,Q, (34) 

+ 00 

K c (x,t,C) = ^(Imro)M(i,t,r ) - ^M(x,t,a ) + % - ^ x(n){M(x,t,T n ) -M(x,t,a n )} , 



2 

n— — oo 



it being understood that all pole or logarithmic singularities, if any, have been removed from M (for brevity, we do 
not introduce special designation for this truncated function). 
To conclude this section, let us show that the velocity field 

$ r +B 

cA*) = ~ J dxK{x,t,g{z)) 

satisfies Q, v (-b + i£) = fl v (b + i£) for all £ downstream, where n v {Q = u v (g~ 1 (Q) d9 ■ Note that since M(x,t,r) 

is real for r real, differentiation of the %-function in the first term in K c yields a purely imaginary quantity canceled 
by its counterpart from K c . This means that when calculating the velocity field, K c depends in effect on z — x + iy, 
while its complex conjugate, on z = x — iy, enabling us to apply the rule (|18p which gives 



d c +B 

n v (() = 9^J b dxK c (x,t,C) 



Therefore, it is sufficient to prove that 

-B r+B 



f-t-a 

dxK c (t,x,-b + i£) = / dxK c (t,x,b + (35) 

-B J-B 



Using Eq. (f2"T|) we write 



X(n) {M(x, t, r n (x, t, -b + - M(x, t, cr n (x, t, -b + %))} 



71— — OO 
+ OC 



= x(n){M(x,t,a n (-x,t,b + it))-M(x,t,T n (-x,t,b + i^))} 

+M(x, t, t (-x, t,b + + M(x, t, n(-5, t,b + 



TL— — OO 
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Similarly, 

M(x, t, to(x, t, —6 + — M(i, i, <To(i, i, —6 + 

= M(x,t,<Ta(-x,t,b + i£)) - M(x,t,Ti(-x,t, b + ig)). 

Adding these identities, changing the integration variable x — > —x on the left of Eq. ([55]) , and taking into account that 
M(— x, t, C) = — M(x, t, C) proves this equality. In fact, it can be shown that fl v (±b + i£) = 0, but we will not need 
this. Thus, the property b + i£) = Q v (b + il;) is a direct consequence of the chain relations (f27)) and antisymmetry 
of M with respect to x. This fact will be used in the next section. 



C. Evaluation of rotational velocity component near the flame front 

Let us now specialize to the case when the observation point z is close to the flame front. In this case, the To-pole 
crosses the real axis of the complex r-plane near its origin. For sufficiently small |r|, one has, according to the 
definition of particle trajectory, 

X(x,t,t — t) = x + w±(x, t)r, 

Y(x,t,t-r) = f(x,t-r) +u + (x,t)r. (36) 
Expanding also fix, t — r) to the first order in r, we find for the complex trajectory 

Z T = z + [w+(x, t) + iu + {x, i)]r, 

so that the equation g(Z T ) = ( becomes 

g(z) + j^{z)[w+{x,t) +iu+(x,t)]T = C (37) 

The value of the rotational velocity component is given by integrating Ko(^,x,t) over x. If x is sufficiently close to 
the point determined by Eq. (f22| . then Eq. ([37]) gives 



T {X,t > O = — ri JZ Tvi . 7=j"W- (38) 



For other x's which are not close to the point satisfying (|22)l . tq has different form depending on the specific structure 
of higher-order terms with respect to r in Eq. (|36p . However, for such x's the pole tq is far from the real axis. On the 
other hand, its exact position in the complex plane is inessential for the purpose of calculating the rotational velocity 
component, because any displacement of this pole adds to velocity terms = 0, as long as the pole does not cross the 
real axis during this displacement. More precisely, it was shown in Sec. lIIIBl that any pole or logarithmic singularities 
are removed from M by the contour deformation, hence, displacements of tq in Eq. (|34p respect conditions a) and 
b). However, changing location of To generally violates c). To ensure that this condition is still met, it is sufficient 
to fulfil the chain relations (|27p. Indeed, we saw at the end of the preceding section that the latter entail equality 
fl v (— b + i£) = Q"(b + i£) for the rotational velocity component. Therefore, the difference of the functions SV(C) 
corresponding to this component before and after the displacement will satisfy c). Thus, we shift the pole To so that 
its new position is given by Eq. (|38[) for all x, and accordingly, shift the infinite sequence of image poles to the new 
positions given by 



/- « C-Abn-q{z) , ^ 4bn + 2b - C - g(z) , , 

j[w + (x,t) + iu + {x,t)\ i\w+{x,t) - iu + (x,t)\ 

Let us check that Eq. ([27ll holds true for the new pole location. Taking into account that —x + if(—x,t) = —z, one 
finds 



T n (-x,t,-b + i£) 
a n (x,t,b + i£) 



—b + it; — Abn — g{—z) 
-9-(-z)[w + (-x, t) + iu + (-x, t)\ 



Abn + b — i£ — g(z) 



j-(z)[w + {x, t) - iu+[x,t)} 
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Differentiating Eq. (p} with respect to x gives 



which together with Eqs. |T]), |2]) prove the first of the two infinite sequences of relations ([27]) . The other is proved 
similarly. Thus, under the assumption of convergence of the series, rotational component of the burnt gas velocity at 
the flame front takes the form 

$ [ +B 

uj v ~ — — Im / dx [x (Im ro) M(x, t, To) — M(x, t, ao) 

2 J-B 

+ x(n){M{x,t, Tn )-M{x,t,a n )}}, (41) 

n— — oo 

it being understood that y is set equal to f(x,t) after the differentiation. 



D. The master equation 

We can finally write the main integro-differential equation relating the on-shcll value of the fresh-gas velocity and 
the flame front position. As we proved in the preceding section, the irrotational component of the burnt gas velocity, 
v p , which was gradually removed from the burnt-gas velocity downstream, satisfies conditions a)-c), so that the 

corresponding function £! P (C) = cj p (g -1 (£)) ~ dQ^ obeys Eq. (fT4f . This relation can be put in a more compact form 

by introducing operator "K defined by 

(fia) (x) = ^|±) j + \xcot{^ h \g{z)-g{z + )])a(x), (42) 

where a(x) is assumed to be continuous function having zero average across the channel, i.e., 

{0) ^ dia{S;) = - (43) 
It is proved in Appendix that under this condition, "K satisfies the usual identity 

"K = -1. (44) 

The dispersion relation for oj p thus becomes 

(l + iA) («$.)' = 0. (45) 

A calculation quite similar to that carried out in Sec. IIII Al shows that the fresh-gas velocity, being irrotational by 
virtue of the Thomson theorem, satisfies 

(l-t3c)(n_)' = 0. (46) 

The proof that the function £!(£) = ^(ff~HC)) ~~^^ corresponding to the fresh-gas velocity fulfils a)-c) (with the 

requirement of boundedness applied upstream) is the same as for Slo given in Sec. IIII Bl Let [v] denote the jump of 
the gas velocity across the flame front, [v] = v(x, f(x, t) + 0) — v(x, f(x, t) — 0). Then the master equation for cj_, / 
is obtained by substituting 

= — — 7+ = 7 (z + ) = — (z + ) , 



Prcmixcd Flame Propagation in Curved Channels, v 1 - October 13, 2009 



12 



in Eq. (|45|). and using Eqs. (|4T1). (146 



7+ 



1 + £K 



7+ 2 7 4 



-Im 



dx 



— B 



X (Im r ) M(x, t, To) - M(x, i, cr ) 



x(n){M(x,t,T„) - M(x, f, cr„)} 



(47) 



By virtue of the identity (|4"4"]) . solutions to this equation automatically satisfy Eq. (|4T)]) . The quantities cr + , u" 
determining the memory kernel, as well as the velocity jumps appearing in this equation, are all known functionals of 
the on-shell fresh gas velocity [H, For instance, in the simplest case of zero-thickness flame fronts one has 



w 



-f- 



N 1 1 ' N 

9-1 ( 2)w_ ,T>u_ 1 T>f 
< h f 1 — 



(48) 
(49) 



where 



T> _ d 
Wt = di 



A 9 



N J dx' 

In addition to that, the local consumption rate in the flame is given by the so-called evolution equation 

(«_,n) = l + S(«_, w-,f), (50) 

where S is a known functional of its arguments, proportional to the cutoff wavelength A c . Together with this equation, 
the complex Eq. (|47[) describes unsteady flame propagation in curved channels in the most general form. 

Let us consider the special case of straight channel. Setting g{z) — z, 7+ = 1, expressions (|39[) can be rewritten as 



T n (x, t, z) — —i-^-(z — [z + 4bn]) , a n (—x, t, z) = —i-^-{z — [z + 4bn + 2b]) . 
v + v + 



where uj + — u + + iw+ , v\_ = w+ + u+ . It is seen that the two infinite sequences of poles span unique sequence 

T n (x, t, z) = — i-^-(z — [z + 2bn\) , n & Z. 

v + 

Furthermore, if we extend the i-domain from [—6,6] to (—00,00), and change x — > —a; in the sum over c-poles, the 
x- integral in Eq. (147[) becomes, taking into account oddness of M and Eq. 



-foe 



dx\ (ImT)M(i,t,T), 



where 



T = T(x, t,z) = —i—^-(z — z). 
Introducing vector r = (x — x, /(x, t) — /(x, t)), this can also be written as 



- - '\::l.nVi-\ 1 



Recalling the definition (jT9j) , the master equation thus reads 





i$ 


f+OC [■ 


){m- 




1 dx 1 












I —00 J r 



dx / drM(x,t-r) 
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where 

T± = 72 ±«— Wl ~ 2 2 • 

This equation coincides with Eq. (5.5) of Ref. [3] where it was shown, in particular, that complemented by the 
evolution equation, it contains all theoretical results on flame dynamics established so far, including the linear equation 
describing Darrieus-Landau instability of planar flames, and the nonlinear Sivashinsky-Clavin equation [14j for flames 
with weak gas expansion. 

To ensure convergence of the improper integral over x, it was regularized in Ref. Q by introducing the factor e~ M '' 
with sufficiently large positive fi, the rotational velocity component being understood as the analytic continuation of 
the regularized expression to fi — 0. A similar but more physical regularization needed generally to ensure convergence 
of the series in Eq. (|4"T)) will be introduced in the next section. 



IV. STEADY FLAME IN A CURVED CHANNEL 



A case of special interest for application of the developed theory is the steady flame propagation. To consider it in 
a curved channel, we have to assume that there is an incoming flow of fresh gas, large enough to hold the front at a 
given position in the channel. Take, for instance, a channel shaped bottlewise as shown on Fig. [3j The map from the 
auxiliary to the physical plane is given by 



- 1 (C) = C + -(B„ 

7T 



6)lnU + exp -— 1 



~2b 



(51) 



where b and are the channel half- width in the far up- and downstream, respectively. If the incoming flow velocity 
exceeds the relative curved flame speed in the straight channel of width b, and B^ is large enough (such that the flow 
velocity at z = +oo is less than the normal flame speed), then the flame will be stabilized somewhere in the convergent 
part of the bottleneck. It is easy to see that the master equation ([47|) still applies to the case under consideration 
despite different initial conditions. In fact, it is evident from the derivation of the expression (141 [) for the rotational 
velocity component that it has the same form whatever the net gas flow. On the other hand, the dispersion relations 
for the irrotational velocity components involve differentiated functions Q = u)/ (dg/dz) which are nothing but the flow 
velocities in the auxiliary plane. Therefore, a constant shift of the velocities does not change these relations either. 
In the steady case, the memory kernel is time-independent, M(x,t) = M(x), so that Eq. (|47|) would take the form 



21^ 

7+ 



+ B 



— — / dxM (x) 

7+ 2 7+ J_ B 



|Im t | - Im ctq + Im ^ ( T « - a n) 



n— — oo 



0. 



(52) 



However, as written this equation has only a formal meaning, because the infinite series in the braces does not converge, 
as can be seen using the explicit formulas (|39[) . One might think that this complication is not essential, since each term 
in this sum is a linear function of £ = g(z). Indeed, if the term-by-term differentiation of the series were allowed, then 
<pg(z) — 2dg/dz taken on-shell would just cancel 7+ in the denominator, leaving us with an cc-independent expression 
which falls off from the equation upon the subsequent x-differentiation. But this operation is not allowed, because 
Eq. (|47J) was obtained under the assumption of convergence of the series. A physical way to resolve this ambiguity is 
to recall that the gas viscosity, though usually negligible in practice, never vanishes in principle. This fact suggests 
the following natural regularization of the series. In the case of a steady viscous flow, Eq. Q is replaced by 

(v, V)er = vAcr, 

where v is the kinematic viscosity. It follows from this equation that for sufficiently small v, vorticity decays along 
the stream lines roughly as a ~ exp(— vr/b 2 ). Since M(x) is proportional to a, this means that the definition (flU|) of 
the function M(i , t, r) is to be replaced by 

+00 



M(x, t, t) = lim / d,TiM(x,t — Ti)e~ 

e^O J 
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(according to the discussion at the end of Sec. HH this integral converges for t\ — > oo, so the limit does exist). The 
parameter e > here is generally a function of x, but for brevity, we suppress this dependence. We need the imaginary 
part of M, which can be written as 

ImM(i,t,r) = — lim / dnM(x, t - Ti)e~ £Tl . (53) 

2l e^O J 

T 

In the steady case, this becomes 

ImM i, t = — ^ lim , 

implying that the series in Eq. (]52fl is to be understood as 

, v.. . , . "I- 

2i e^o e ' L J 




, summing the geometric series, and taking the limit e — > gives 



[C-.g(5)] 2 , [C + .g(5)] 2 



D J 

where D = -D(i) = ; y[w+(x) + iu+(x)\. It should be mentioned that the limits of the sums over r- and cr-poles do not 
exist separately. For instance, one has 



+oo , s 2 cosh 

X(n)e~^ =exp|j[C + ^l)]| 



r2e& 



^ -=- !> 1 



(54) 



It is evident that the imaginary part of this expression is 0(1 je). It is canceled by a similar contribution coming from 
the sum over T-poles, as is the next 0(l)-term. In effect, the exponents in the numerator are to be expanded to the 
second order in s, which is the reason why the resulting expression turns out to be quadratic in £. Putting it into 
Eq. ([52]) . and performing ^-differentiation with the help of the rule l|T8)) yields 



•>{ U ~X j-fij-'v}) M P B **Tur(-\ r x+(Imr ) g(z+) ( 1 1 V " ' 

Finally, noting that x+ (I mT o) — ~ x )i the x-derivative of the first term in the integrand removes the x-integral, 
while changing x — > — x in the second and taking into account that D{— x) = —D(x) brings the master equation to 
the form 



= 0. 



-V + (i + &)< ( M V - M(a) r; ( r } + ^¥ r ds*?*;® 

l+J v J \\1+J l+v+(x) 2b J_ B 7<W 



The last term in this equation, resulting from the infinite series over r- and cr-poles, ensures that the expression in 
the braces has zero average across the channel, as it should be according to the definition of the operator J£ . Indeed, 
it follows from g(z+)\ x =±B = ±b + i£ that the average of the last two terms in the braces vanishes, while 

1 [w\ 




2B 74 

Evidently, the quantity I/7+ taken at the channel wall represents vector normal to the wall at the front end-points. 
Since at these points gas velocity is tangent to the wall, so is its jump, which means that [w]/7 + is real there. On the 
other hand, the values of [w]/7+ at the left and the right wall are complex conjugates of each other, as is seen from 
Eqs. ©, (|4"D|) , thus proving that the average of the first term in the braces is also zero. 
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2(u J .y + (l + m){[.}'-^^ + f± I ^h0, (56) 



Let us consider flame propagation in a straight channel in more detail. The master equation reduces in this case to 

M(x)w+(x) z' + r+ b M(x)w+(x) 
v\{x) + 2b]_ b dX v\{x) 

This equation differs from the corresponding equation of Ref. [5j exactly by the last term. The point is that in 
Ref. [5J, channel flame propagation was considered as the propagation of a 26-periodic infinite flame; though not 
stated explicitly in Ref. vanishing of (Muj + /v\) is actually the condition of existence of improper integrals over 
the infinite flame front, appearing in the derivation of the rotational velocity component. The question whether or 
not this quantity is zero turns out to be quite nontrivial. As to its real part, it is zero identically as the average 
of an odd periodic function. However, there is no apparent reason why {Mw+/v\) should vanish. Using the small 
{9 — l)-expansion, it can be shown that this quantity does vanish at least to the fourth post-Sivashinsky approximation 
(which corresponds to keeping 0((9 — l) 7 )-terms in the master equation), as an average of the derivative of a periodic 
function [l5j . Furthermore, it vanishes in the opposite case of flames with large front-slope (such as flames anchored 
in high- velocity streams), which can be shown even without using the master equation, but we were unable to prove 
this in the general case. It is interesting to note that the new term, if nonzero, would mean appearance of an effective 
gravitational field. Indeed, in the simplest case of constant gravity parallel to the channel walls, a term proportional 
to /' appears in the on-shell vorticity [lj|, which has the form 

Aa + = -^p-gf'(x,t), 

where g is the gravity acceleration. Since z' + — 1 + if'(x, i), it follows from Eq. (|5"6"|) that taking roughly u + ~ v + ~ 6, 
the effective acceleration can be estimated as g e ff — —(Mw+/ig_). In this connection, it is worth of mentioning that 
the extra term belongs to the kernel of the !K-operator, namely, in the general case of curved channel, "Kg'(z + ) = 0, 
as can be easily checked using the definition 



Quasi-steady flames 

Characteristic time of spontaneous evolution of a flame disturbance with characteristic wavelength A is, by the 
order of magnitude, T ~ A/t/f, where Uf is the flame front velocity relative to the fresh gas. Hence, dynamics of 
disturbances with sufficiently large A is slow, and the flame can be considered quasi-steady, provided that the steady 
regime exists in the absence of disturbances. Specifically, since there is only one internal time scale for steady flames, 
namely, A c /C/f , such consideration is applicable, if A 3> A c . Flames propagating in channels with slowly varying width, 
or in a slowly varying gravity deliver simplest examples of quasi-steady flames. The results obtained above for steady 
flames arc readily generalized to this case. 

In terms of the Fourier decomposition of the memory kernel, 



M(x,t) = / diuM{£,uj)e- luJt , 

J — oo 

slowness of flame dynamics means that M(x,oj) is noticeable only for sufficiently small \oj\. Substituting this into 
Eq. (I53p , the e-regularized expression for Im M can be written as 



1 f + °° f 

ImM(x,i,r) = — lira / duM(x,uj)e- lut / dne^-^ 1 

2lz->°J-oo J 



1 r+°° p -(e-iu))r _ p —(e-iu)T 

= — lira/ du)M(x,cj)e- %ut ; . (57) 

It is seen that as a result of the flame unsteadiness, the regularizing parameter e becomes effectively complex. To 
sum the series in Eq. (|47[) in the steady case, we had to expand the right hand side of Eq. (|54[) in powers of e keeping 
terms 0(e). Therefore, in order to take into account the unsteady effects to the first order in uj, we have to retain in 
Eq. (|5"4"|) also the term 0(e 2 ), and to replace e — > — iui afterwards. A simple calculation gives 

n— — oo 
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where now D = D(x, t) = 7 [10+ (a;, t) + iu+(x, t)], and contributions linear in £ have been omitted (as before, they fall 
off from Eq. (|T7|) anyway). Hence, 



+00 



Im x(n){M(x,t,T n )-M(x,t,a n )} 



n— — 00 



(~iu) 



126 




1 dM(x,t) 1 
~ 2i 9* 126 

With the same accuracy, one has 



{...} 



lm.M(x,t,T ) = — 

2,1 



M(x,t)(n-T ) + 



aM(i,t)r 2 -7s 2 

at 2 



and a similar expression for Im M(i, t, uq). Substituting these formulas in Eq. (|47|) . and performing differentiations as 
before yields the master equation for quasi-steady flames 



7+ 



1 + 




M(x,t)w+(x,t) , g'(z+) f +B M(x,t)w+(x,t) 



7+ 



2 



B d£ / u+(x,t) y dM(x,t) 

b x \iv\{x,t)j at 



21, J_g W+(x,t) 

g( z +) - g{z) 



X(x -x) + 



= 0. 



It is straightforward to check that the new term has zero average across the channel, as it should. 

V. THE SMALL GAS EXPANSION LIMIT 

A. Derivation of equation for the flame front position 

Equation (|4T[) also considerably simplifies in the case when (9 — 1) = a is small. Since the cutoff wavelength 
A c = 0(1/ a), the wave numbers of the flame perturbations are 0(a), implying that spatial differentiation of a 
perturbed flow variable raises its order with respect to a by one. Next, estimation of the perturbation growth 
rate using the Darrieus-Landau relation shows that it is 0(a 2 ), therefore, time-differentiation adds two orders. The 
evolution equation ([50]) [with S = 0] then tells us that u_ = 0(a 2 ), hence, W- — 0(a 2 ). Also, the slowness of flame 
dynamics means that the channel width must vary sufficiently slowly along the channel, namely, the phase of the 
function 7(2) = dg(z)/dz must be small, arg7(z) <C 1. A natural representations of the maps meeting this requirement 
is 



// : {j) 



(58) 



where L = 0(b/a 2 ) 3> b is the characteristic length scale over which the channel width changes significantly, and the 
function h(£) is "normal," in the sense that \dh/dC\ = 0(1) [dh/d( denotes hereon derivative of the function h(Q, not 
h((/L)]. It obeys the usual reflection rule h (— C) — — h (0 • It follows from this relation that the derivative of g~ 1 (C) 
given by Eq. (|58[) is a real 0(l)-quantity at the imaginary axis (£ = which acquires an 0(b/L) imaginary part at 
the channel walls (£ = ±6 + Using this representation, the factor 77/ appearing in Eq. (jlT|) can be evaluated as 



1 

7+ 



dq 1 . „ . dh , „ / t \ dh , _ d 2 q 1 , . . . > 

-^(C) = Tc (C + /l) = Tc WL) + -|^o)(C, 



»&) + O 



(59) 



where £0 1S the image of the middle-point of the flame front, i£o = <7(i/(0, t)). The quantities 

dh , , . , . 1 d 2 q~ x , 
—(i^/L) = s(t), ? -^5-(i£o) = k(t) 
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are thus purely real. We note for future reference that they can be related to the flame front velocity at x = 0, namely, 
differentiating s(i), and using the above definitions gives 

ds(t) d 2 g-\ dg , k(t) df(0,t) k(t) 4 

The last expression follows from the fact that the flame propagates downwards, its velocity being 1 + 0(a 2 ). 

Below, we consider the first post-Sivashinsky approximation which corresponds to retaining terms 0(a 4 ) in Eq. f|47[) . 
Taking into account that cr+ = 0(a 4 ) [Cf. Eq. (05])], one has to this order 

M = a + = -a(?£- + f'f 



, 9t 

Moreover, since M is already of the fourth order, r„, er n can be taken in the form 



T n = is(t)[4bn - C + 9{z% °n = is{t)[Abn + 2b~Q~ g{z)}. 
One has, within the given accuracy, 

M(x,t,r n ) - M(x,t,a n ) = / dn(7 + (x,t - n) = a + (x,t)(a n - t„) 



This is independent of z, which means that the series in Eq. (|47[) requires no care, and the image contributions can 
be omitted altogether. On the other hand, for the ro-contribution one finds ImM(i,i,To) = t)s(t)[£ — g(z)]. It 

is worth of mentioning that the possibility to replace M{x, t — r) by M(x, t) signifies the absence of memory effects in 
the present case. When expanding various terms in Eq. (|47p in powers of a, it is sufficient to retain only the (9(l)-part 
in the expression (|59p . the only exception being the contribution proportional to [u]. Namely, using Eq. (|48p one has 



JUj 

7+ 



= ^ [s(t) + ik(t)(C+ - i^o) + 0(a 4 )] - as(t) (l - ^\ + iak(t)r, + 0(a 4 ), 



where we took into account that — = 0(a). Writing also r/ = x/s(t) + 0(a 2 ), substituting this into Eq. (I47p . 
and carrying out differentiations as in the preceding section gives 



2(.->< + a( 1+1 *){«J§-,f + f}' 



0. (61) 



Next, it follows from Eq. (|59p that the value of I/7+ at the right channel wall is s(t)[l + ik(t)B/s 2 (t)]. Considered 
as a vector in the complex z-plane, this is normal to the wall, hence, xk(t)/s 2 (t) — f vanishes there. Therefore, the 
expression in the braces in Eq. (|6ip has the same value for x = ±B. For such a function a(x), the !H -operator acting 
on its derivative can be rewritten as 



(%a'){x) = 9 -^j b dx cot {^[g(z)-g(z + )}} a' (x) 

B ( n 1 

dxg'(z) cot { — [g(z) - g(z+)]j a{x). 



1 d 
2bdx 



-B 

Furthermore, g'(z) = (dg/dz)z' can be replaced here by z'/s(t), with an error of the order 0(a 2 ). Noting also that 
within the same accuracy, B — s(t)b, we find 

{x) = hi tl d2cot {jb [ "- z+] } a(i) - (62) 

But the latter expression can be shown to be just (H B a)', again up to terms of the second relative order in a, where 
Hb is the usual Hilbert operator defined on the interval (—B,+B) (Cf. Eq. (B12) in Ref. ||) 

(h b oJ (x) = i j- dxcot\J^[x -x]}a(x). (63) 
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Integration of Eq. (|6"T]) thus yields 

2._ +Q (l + iJ J fl ){ i ^- i / + f}= C(t ), 

where C(t) is an "integration constant." The parity properties © imply that it is real. Taking into account that the 
Hilbert operator is also real, extraction of the real and imaginary parts of this equation gives 

2u. + aH B ^f-x^\+Jl = C{t), (64) 
2w_ - af = 0(a 3 ). (65) 

The latter equation is written out up to 0(a 3 )-ter ms, because the ui-component of the fresh gas velocity is multiplied 
by /' = 0(a) in the evolution equation (f50|) which in the present case reads 

df f' 2 

«- - fw- = ^ + 1 + V • ( 66 ) 

Combining the last three equations, we arrive at the following equation for the flame front position 

(0 + 1)§[ + 9f 2 = -(0 1)H B {/' - xMj + C (t) 2. (67) 
It is convenient to introduce a new function 

F(x,t) = f(x,t)-'^, F'(±B,t)=0, 
in terms of which Eq. (|6T[) takes the form, within the third order accuracy, 



1)^ + OF' 2 + = -(» - I > //« F' O / 1 - 2. ( (iM 



Averaging this equation along the front yields the function C(t) itself 

C (t) = 2 + (9 + l)^- + 6(F' 2 ) + ^(xF'). 

In the case of a straight channel, k(t) = 0, and the obtained equation is simply the Sivashinsky-Clavin equation 
[l4| (with the term C(t) added according to 17]). Equation (|6^|) is valid for zero-thickness flames. Account of the 
transport processes inside the flame front adds higher-derivative terms proportional to the cutoff A c . If, as usual, these 
processes can be considered linear, then the term to be added to the right hand side of Eq. (|6"?| has the form [l(| EH 
(6 — l)A c /"/27r, and hence, Eq. (|6"8|) is to be replaced by 



(0 + 1) 

where 



dF ds xF' 
~dt + dts^j 



OF' 2 = -{0 - 1)H B F' + {6- 1)^F" + C(t) , (69) 



2tt s 2 (t) ' 

and we also used Eq. ([60]) to rearrange the left hand side. 

Concerning the structure of the obtained equation, it is important to note that although the function s(t) is related 
by its definition to the front location within the channel, it is actually independent of the specifics of flame evolution. 
Namely, it can be written as 

s(t) = § ( ^fo -^) ) t {70) 



d( V L 

where fo = /(0,0). Since the flame velocity increase due to the front curvature is 0(a 2 ), corrections to this formula 
are of the fourth order. Thus, in the first post-Sivashinsky approximation, the function s(t) is to be considered as 
non-dynamical quantity predetermined by the channel geometry and initial front position. 
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B. Front evolution in terms of pole dynamics 



Though not obvious from its present form, Eq. (|69|) is amenable to the pole decomposition. To see this, let us 
introduce new variables 77, $ according to 



F(x,t) = $ 



e(t) 



a(t) 



Then the expression in the square brackets on the left of Eq. (f6"9l is just d^/dt, so that this equation becomes 



0£ 

1 aj 



dt s 2 {t) \dr) 



<9$ 



-'■4? 



- 1 A r d 2 $ 



s(t) dr] s 2 (t) 2vr drf 



+ C{t) , 



(71) 



where Hi, is defined by 



(H b oj (77) = ^ j dr) cot I ^[V-V] } a iv)- 



Solutions of Eq. (J7TJ) can be found in the form [l8|, [l^| 



2P 

$(77, t) = * (t) + A £ In sin { ^ - % (t)] } . 



(72) 



fe=i 



Here P is the number of complex-conjugate pairs of poles located at the points r)f-, k — 1,...,2P, in the complex 
?7-plane. To determine the amplitude A and the functions <&o(t), Vk(t), we substitute this anzats into Eq. (|71j) . and 
use the formulas 



ip 

^ = -^E^ + ^Imr/fe) cot {^[77-77^*)]}), 



k=l 



cotscoty = — 1 + cot(x — y)(cot y — cot x) 
This leads to the following system of ordinary differential equations 

- 1 A r 



A = — 



Mo 
dt 

dr/ k 
dt 



9 2ir ' 

- I) 2 PA C / PA C 



' + 1) 2bs(t) \2bs(t) 
- 1 1 



1 



C(t) 

(e + iy 



ls(t) 



IP 

^ Im%) + 2^)g C ° t {| (?7fc " 77ro) } 



(73) 

(74) 
(75) 

(76) 



As is seen, this system does not determine the functions 3>o(t)> C(t) individually, but only a linear combination thereof. 
This corresponds to the fact already mentioned above that the master equation has one and the same form whatever 
velocity of the incoming fresh gas (potential) flow. We now recall that the fresh gas is assumed to be initially at rest. 
The flow continuity then implies that its total flux through the front is zero. Using the evolution equation (|66p . this 
condition can be written as 



dt 



</' 2 > 



In terms of the original variables, Eq. I|72p reads 

x 2 ds 



2s(t) dt 



ip , 
+ A£lnsin||- 

fc=i k 



-77T -Vk(t) 

s(t) 



(77) 



(78) 
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Substituting this into Eq. (|TT[) gives, within the given approximation, 



M Oh Hi O R 



k=l 



2b 



1 / ttA 



1 



2 \2bs 2B 



2P 



-B 



cot 



\/c=l 



a; 





2/* 



(notice cancelation of terms proportional to ds/dt). A straightforward calculation using the formulas ([75]) . (|7d| . and 



gives 



eft 



= -1 - 



>-l) 2 PA C 
26 2 2bs(t) 



PX C 



[ - I) 2 A c 
26» 2 26 



E 



rflm ?? fc 
dt 



where the latter sum runs only over the poles in the upper half-plane. Integrating this relation with the same accuracy, 
we find 



$o(t) 



-t - 1- 



>~1) 2 P\c 

26 2 2bs{t) 



2h,,(t)) 



I) 2 A c 



29 2 2b 



Im rjk (t) + const, 



fc=i 



where the additive constant determines initial front position. 

To summarize, in order to determine the flame front evolution in a channel, one uses the given map g _1 (C) and the 
quantity /o to construct the function s(t) according to Eq. (|70p . and then solves the ordinary differential equations 
(|76f for rjk(t). Substituting the latter into Eq. (|78|) yields the function f(x,t). An important point of this procedure 
is the proper choice of the number P. Since by the construction the number of pole pairs is independent of time, any 
choice of P predetermines to a certain extent the flame evolution. At the same time, it is known that in straight 
channels, there is only one stable steady pole solution for every channel width [2fJ, [21[. The corresponding number of 
pole pairs is given by 

( B 1 
Ps=Int U + 2 

where Int(x) denotes the integer part of x. Thus, in order to construct the true solution which correctly describes 
unstable flame behavior in curved channels, P is to be taken large enough. In channels of finite width, it is sufficient 
to set P s = Int (-B m ax/Ac + 1/2) , where 2_B max is the maximal channel width. Then in regions with B < -B max , a 
number of redundant poles move to infinity in the ?7-plane. In diverging channels, however, P s is formally infinite, so 
that this case requires special consideration. 



VI. DISCUSSION AND CONCLUSIONS 

The results obtained in this paper extend the on-shell flame description constructed in Refs. [1, S 13] to the case 
of flames propagating in curved symmetric channels. Namely, complex Eq. (|4T|) relates the on-shell gas velocity and 
front position of thin unsteady flame in the most general form. This equation shows that the influence of channel 
geometry on flame dynamics is deeper and more diverse than might be expected from simple considerations based on 
potential flow models [10|, • In the latter, the effect of wall curvature is contained entirely in the quantity 1 /j + 
accompanying w_ and [ui] in Eq. ([TT]) . which is just a conversion factor relating on-shell velocities of the potential flows 
in the physical channel and its image. The master equation shows that things are actually much more complicated. 
In addition to the usual nonlocality related to potential flows and described by the operator Jf , vorticity produced 
by the curved flame brings in specific nonlocality associated with the integral term in Eq. (|Tf)) . The structure of this 
term essentially depends on the properties of the map g(z). Moreover, temporal nonlocality described by the memory 
kernel is also affected by these properties through the r- and c-poles appearing in the function M(x,i, r). 

It is worth of recalling that the derivation of Eq. (|4T|) relies on two important technical assumptions, namely, 
convergence of the integrals over r corresponding to the remote past, and convergence of the series over the r- and 
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tr-poles in Eq. (j34|) . Concerning the former, asymptotic straightness of the channel is a sufficient condition, as was 
discussed in Sec. [D] but more general channel configurations may require additional consideration such as for instance 
stability analysis of a particular jet configuration. As to the series convergence, it seems to pose no serious problem 
as it can be naturally ensured on physical grounds. Indeed, as was demonstrated in Sec. IIV1 even in the case of a 
steady flame when assumptions about the large-time behavior of the memory kernel are irrelevant, convergence is 
easily achieved taking into account small gas viscosity. 

Finally, we mention that flame anchoring systems such as metallic rods can be easily incorporated into the above 
construction following the lines of Ref. [8j . In particular, this gives immediate analytic access to the problems involving 
fast flow burning, because the IK-operator considerably simplifies in this case (see Refs. @, 0] for details). 
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APPENDIX. PROOF OF THE IDENTITY "H 



Let us prove that "K satisfies the operator identity "K = — 1 valid on the space of functions a(x) which are analytic 
in a; in a vicinity of the real axis, and have zero average across the channel (the proof follows closely that given in 
Ref. [j| in the case of straight channel). For this purpose we rewrite its definition as 



"Ka) (x) 



g'(*+) 
46 



dxi 



cot {^[ff^i) ~ 9(z+)]j a ( x i), z i = x i + i f( x i, t ), 



(79) 



where C\ = C a U is the contour in the complex ir-plane, shown on Fig. [4] We then have 

(?{ 2 a) (as) = (?C(5fa)) (x) 
<?'(*+) 



46 



] - J dxi cot {^Wi) - 9(z+)}} J dx 2 cot {^{g(z2) - g(zi)}} 

C\ C2 



a(x 2 ) , 



where z 2 = x 2 + if(x 2 ,t), and the contour C 2 = U of integration over x 2 comprises C\ (see Fig. |4j. The 
contours Ci, C 2 are chosen so that all singularities of the integrand (if any), except the poles of the cotangent, remain 
above C^ 1 ", or below C^~, which is possible under our assumptions about analytic properties of a(x), g(z). Now we 
change the order of integration, and use the formula (|73|) to compute the integral over x\ 

dg{zi) 
46 



cot {^[ ff<Zl ) -3(^+)]}cot{^[ 3 (z 2 ) -g(zi)]] 



Ci 



din- 



sin 



{^ h Wi)-g&)]} 
{^Wi)-9{*+)\} 



In view of the 26-periodicity of the integrand, it is phase change of the argument of the logarithm that only contributes. 
It is equal to 2iri when x\ runs and x 2 € ■ (—2iri) when x\ runs and x 2 S C 2 , and zero in the other two 
cases. Taking into account also Eq. (|43|) . we thus find 



(lK 2 a) (33) 



ig'(z+) 
46 

ig'{z+) 

46 



ig'{z+) 
46 



dx 



10(12) cot {^[9(22) -g( z +)}} 



dx- 



>a(x 2 )cot|^[3(z 2 ) -g(z + )]j 



2m ■ res 



a(x 2 ) cot \ ^-{g(z 2 ) - g(z + )]\ = -a(x) 

I ZD J xo—oo 



(80) 
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as was to be proved. 
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